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D , ABSTRACT 

In the nonrelativistic case I found that whenever the relation mc 2 /e 2 < X(a,g m ) is 
satisfied, where a is a flux, g m is magnetic moment, and X(a,g m ) is some function that 
is nonzero only for g m > 2 (note that g m = 2.00232 for the electron), then the matter is 
unstable against formation of the flux a. The result persists down to g m = 2 provided the 
Aharonov-Bohm potential is supplemented with a short range attractive potential. I also 
show that whenever a bound state is present in the spectrum it is always accompanied by 
a resonance with the energy proportional to the absolute value of the binding energy. In 
the relativistic case one again finds the resonance when the bound state is present but the 
instability disappear as long as the minimal coupling is considered. For the Klein-Gordon 
equation with the Pauli coupling which exists in (2+1) dimensions without any reference 
to a spin the matter is again unstable for g m > 2. The results are obtained by calculating 
the change of the density of states induced by the Aharonov-Bohm potential. The Krein- 
Friedel formula for this long-ranged potential is shown to be valid when supplemented 
with zeta function regularization. 
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In this letter a change of the density of states in the whole space is calculated for the 
Schrodinger, the Klein-Gordon, and the Dirac equations with the Aharonov-Bohm (AB) 
potential A r = 0, = $/27rr. This enables to discuss the stability properties of matter 
against the spontaneous creation of a magnetic field. In the relativistic case and for 
the normal magnetic moment we reconfirm previous result, known as the diamagnetic 
inequality 0J that the matter is stable. The latter was proven under the assumption of 
minimal coupling which implicitely assumes the normal magnetic moment. However, in 
the nonrelativistic case a window exists for the magnetic moment g m > 2 in which the 
inequality is violated leading to the instability of matter against a magnetic field formation. 
The reason is the formation of bound states which decouple from the Hilbert space by 
taking away negative energy. 

Note that one has the unitary equivalence between a spin 1/2 charged particle in a 2D 
magnetic field and a spin 1/2 neutral particle with an anomalous magnetic moment in a 
2D electric field ||. One also has a formal similarity between the scattering of electron in 
the AB potential and in the spacetime of a gravitational vortex in 2 + 1 dimensions M. 
Moreover, one encounters the AB potential (of non-magnetic origin) in the cosmic-string 
scenarious || and our results apply to the this cases as well. 

1. The Schrodinger equation.- By using the separation of variables, assuming e = — |e|, 
the total Hamiltonian is written as a direct sum of channel radial Hamiltonians Hi || ||| 

s 1 d 1,2 
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where v = \l + a\, a being the total flux $ in the units of the flux quantum $ Q = hc/\e\, 
with the spectrum given by 

Mr,<P) = J\i +a \(kr)e^, (2) 



with k = ^j2mE/h 2 . The eigenfunction expansion for the Green function in the polar 
coordinates x = (r x ,(p x ) is 



G(x,y,£) = -W £ e^-^JiHaiMJit+aiCAr,). (3) 
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The normalization is chosen such that — (l/7r)ImTr G(x, x, E + ie) gives the two dimen- 
sional free density of states p Q (E) = (mV/27rh 2 )dE for a G Z. In this case the sum in 



1 



(0) can be taken exactly by means of Graf's addition theorem (||7), relation 9.1.79). By 
taking the integral (assuming q 2 = q 2 + is, and r x < r y ) one finds 

G(x - y - E) = w L~ j^'-M* - = <«i* ~ y|) - (4) 

One can also take the integral 

roc hdk jci 

/ TT TT 1 M kr x)Ju{kr y ) = - — J u (qr x )Hl(qr y ) (5) 

Jo q -\- %t — K a 

at first (which is valid for nonintegral v too) and then to take the remaining sum by Graf's 

theorem with the same result (Q). From the 'formal scattering' point of view taking the 

residuum at k = \q\ + ie corresponds to choosing the outgoing boundary conditions. 

From the point of view of L 2 (R 2 ) it corresponds to taking the boundary value of the 

resolvent operator on the upper side of the cut at [0, oo) in the complex energy plane ||. 

The limiting value of the resolvent operator on the lower side of the cut is the complex 

conjugate of (|J), the discontinuity across the cut G(x, x, E + ) — G(x, x, EJ) = —i(m/h 2 ), 

and 

- (l/7r)ImG(x,x,£ + ) = m/2iTh 2 (6) 

which confirms our normalization. 

Whenever a $ Z Graf's theorem cannot be used. To proceed further with this case we 
use the fact that @ has an analytic continuation on the imaginary axis in the complex 
momentum plane 

g2 + k2 Jv{kr x ) J v {kr y ) = —J u (iqr x )Hl(iqr y ) = I^qr^K^qry), (7) 

where I u and K v are modified Bessel functions. To sum over I one uses the integral 
representation of these functions and following the steps of f| one can separate the in- 
dependent contribution 

msinM /•« r°° , .^Mr.rcoshtf+coshu,) e^-") , o , 



G Q (x,x,M) -G (x,x,M) = r M f° ^-^.(c-h*^-) 

Ti J-od J-oo 1 + e 
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where M = —iJlmEjU and r/ is the nonintegral part of a, < rj < 1. After taking the 



trace over spatial coordinates, using formulae 3.512.1 and 8.334.3 of [HJ, and returning 
back on the real momentum axis one finally finds 

Tr [G a (x, x, E) - G Q (x, x, E)] = ~ V (l - V ) ± (9) 



which gives the change of the density of states in the whole space 



p a (E)-p (E) = -^rj(l-rj)6(E) 



(10) 



that confirms JTTJ] , where it was obtained in the context of anyonic physics. 

The channel Hamiltonians Hi for which 1 1 + a | < 1 admit a one-parametric family of 
self-adjoint extensions P, RJ. We shall consider the situation with bound states 



(11) 



of energy E\ = — (h 2 /2m)K 2 at I = —n and I = —n — 1 channels, with n = [a] denoting 
the nearest integer smaller than a. In the presence of the bound states the scattering 
states (@) have to be modified 



ipi(r,4>) ipi(r, 4>) = (J\i +a \(kr) - AiJ^\i +a \(kr))e % 



(12) 



This is because Hi has to be necessary a symmetric operator what already determines 
Ai to be Ai = (k/Ki) 2 ", i.e., energy dependent. To calculate the change of the integrated 
density of states in the whole space we make use of the Krein-Friedel formula |L2| 



N a {E) - N (E) = -— IndetS, 
Ztt 



(13) 



S being the total on-shell S-matrix. Despite the AB potential is long-ranged we have 
found that the Krein-Friedel formula when combined with ^-function regularization can 
still be used. The radial part of the general solution ([12]) behaves as 



Ri(r) 



( 1 _ A,p in \ l + a \ N 

const e~ lkr + 1 . „, i e -i<H<x\+i/2) e ikr 



1 - Aie- in \ l+a \ 

which determines the l-th channel S/-matrix to be S; = e 2tSl with 



r — ► oo 



5i = 2 7r d^l — 1^ + a \) + arc tan 
Without the presence of bound states (Ai = 0) 

oo 

lndetS= 2iS i 



sin(|Z + a\ir) 



^cos(|Z + a\ir) — A 
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in [2C R {s)-C H (s,r])-CH(s,l-v)}\ ( 
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where (r and (h are the Riemann and the Hurwitz (^-function. The use of the Krein- 
Friedel formula (|13D then independently confirms fllPp. Under the presence of bound 
states the contribution of scattering states to N a for E > is 

N a (E) - N (E) = - -r](l + ~ arctan ' s^W) \ 



'/ 



2' v " 7T \cos(r]TT) - (\E_ n \/E) 

1 / sin(?77r) \ 

" ^ arctan Uo^) + (i^-ii/^) (i ^v' ( J 

where E_ n and £?- n -i are the binding energies in I = —n and / = — n — 1 channels. Note 
that for < rj < 1/2 we have a resonance at 

£ = 7 /""Ai/ > 0- ( 18 ) 
(cos(??7r)) 1 /'? v 7 

The phase shift 8- n (E) (0) changes by 7r in the direction of increasing energy and the 
integrated density of states ( [L7|) has a sharp increase by one. For 1/2 < r] < 1 the 
resonance is shifted to the I = —n — 1 channel, rj = 1/2 is a special point since resonances 
are in both channels at infinity. Since we always are in front of the resonances the 
contribution of arctan terms in ( |T7j ) does not vanish as E — ► oo but instead gives —1. 
Different self-adjoint extensions correspond to different physics inside the flux tube 



| I3| . To identify the physics which underlines them we have considered the situation 
when the AB potential is regularized by a uniform magnetic B field within the radius R. 
In order for a bound state to exist the matching equation for the exterior and interior 
solutions in the Z-th channel, 

^Wj _ , ,,| , |/| + / + 1 + (x 2 /2a) iF\ ggg + (x 2 /4a), \l\ + 2, a) 
X K ll+al (x) a + ^ +a |/| + i ^(^i + ^a),!/! + !,«)' 1 j 

with i-Fi(a, 6, c) the Kummer hypergeometric function 0, has to have a solution x\ = 
KiR 7^ 0. However, since the l.h.s. decreases from — \l + a\ to — oo as x —>■ oo and the 
r.h.s. is always positive one finds that it is impossible unless it is an attractive potential 
V(r) inside the flux tube, 

V(T)\r< R = ~-^c{R), (20) 



with c(R) = 2(1+ e(i?)), e(R) > 0, and e(R) -> as R -> @. This amounts to changing 
x 2 /2a by x 2 /2a — c/2 on the r.h.s. of flHf). The attractive potential can be either put in 
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by hands or, when the Pauli equation is considered, as arising in the magnetic momentum 
coupling of electrons with spin opposite to the direction of magnetic field B. The case of 
gyromagnetic ratio g m = 2 (e(R) = 0) is a critical point at which ([19]) has a solution at 
x = for I = —n. It is known that there are ]a[— 1 zero modes in this case, }a[ the nearest 



integer larger than a ||15|| . The result does only depend on the total flux a and not on a 
particular distribution of a magnetic field B. 

Whenever g m > 2 (or g m = 2 with an attractive potential V(r) = —e/R 2 , e > 
arbitrary small) the bound states may occur in the spectrum. They correspond to solutions 
Xi > of ([11]). In contrast to the zero modes their number does depend on a particular 



distribution of the magnetic field B but it is less than or equals to 

# 6 = l+n + 2[a((/ m -2)/4] (21) 

with [.] as above. The bound is saturated when one uses the cylindrical shell regularization 
of the AB potential |T^, [T7|. Note that in this case the energy Eb of magnetic field is 



infinite for any R ^ in contrast to homogeneous field regularization when 

1 <f> 2 

E B = -ttB 2 R 2 = — — • (22) 

2 2tiR 2 v ; 

2. Energy calculations - Up to g m = 2 no bound state is present in the spectrum and 
the change of the density of the scattering states is still given by fllQ|) . Zero modes which 
may occur for g m = 2 are regular at the origin and do not change phase shifts except 
for I = — n and I = — n — 1 channels where they cause phase shift flip (Af 1 = in (|15ft) 
(cf. []T7|]) but in sum they cancel. Therefore the energy of magnetic field (^2]) tends to 
infinity as R — > (when the total flux $ is kept fixed) the matter is stable with regard to 
a spontaneous creation of the AB field. 

When g m > 2 then bound states may occur in the spectrum. Their energy is 

e ' - < 23 > 



and tends to — oo as 1/R 2 when R — > 0, in the same way as the magnetic field energy ( f22"D 
goes to oo. The bound states decouple in the R — > limit from the Hilbert space L 2 (R 2 ) 
(Ai — > in (|1"2D in the limit) and take away the nonperiodicity of the spectrum with 
regard to a — > a ± 1 which persists for any finite R. What is left behind is nothing but 
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the conventional AB problem with the change of the density of states (10). We recall that 
to have a finite energy bound state in the spectrum it has to be the attractive potential 
(POD inside the flux tube. 

Bound states solutions x\ for the homogeneous field regularization determine the func- 
tion X(a, g m ), 

X(a, g m ) = £ x 2 (a, g m ) > 0. (24) 

By comparing the coefficients in front of 1/R 2 in ( p2|) and (|23|) one finds that whenever 

mc /e 2 < X( (25) 

the total energy of field and matter altogether goes to — oo as R — ► 0. 

3. T/ie Dirac and the Klein-Gordon equations - A connection with the Schrodinger 
equation is established by noticing that for the Dirac spinor to belong to the spectrum its 
up/down radial component has to be an eigenfunction with an eigenvalue k 2 of 

d 2 1 d v 2 a „, . . . 

H l = -—--- + -± g m -5(r), 26 
dr z r dr r z r 



with g m = 2 and (z/ = l+a/u = l+l+a) for the (up/down) component ||17fl . The dispersion 
relation for k is now changed to k — (l/hc)y/ ' E 2 — m 2 c 4 . The two-component solutions 
of the massive Dirac equation have only one degree of freedom which is reflected in the 
equality of up and down phase shifts J|. For the conventional AB problem 5f = 5f = na/2 
if I = —n — 1, and 8f = 5f = n(\l\ — \l + a\)/2 otherwise. Using the results of JTJ, |18 



one 



finds that for a fixed real k ^ the spectrum is symmetric with respect to the origin, in 
accord with supersymmetry |TH]. Thus the contribution of scattering states to the density 
of states is given by 

p a (£) - Po {£) = ~ V (1 - V ) 5(£ 2 - m 2 c 4 ). (27) 



The only asymmetry can occur for the threshold [0 and bound states. Depending on the 
orientation of B the threshold states occur either only at the upper or only at the lower 
threshold [|18|]. Their number equals to the number ]oc[— 1 of zero modes of the Schrodinger 



equation 15, 18 



By considering different self-adjoint extensions one finds that, in contrast to the 
Schrodinger equation, bound states can only occur in a single channel I = —n — 1 (I = —n 
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for positive charge) fl3] . In the presence of the bound state of energy E in the I = —(n+1)- 
th channel, 

»*-^,(t,r,#) = i f y^±|^M ) e -.(»«)*e-'"/», (28) 
iV \ iymc z — EK v {Kr)e v J 

where k = (l/h^y/m 2 ^ — E 2 , the scattering state with energy £ > mc 2 is 

^.-n-xft r, 0) = ( ) g-i^iWe-^/ft, ( 29 ) 

v (r\ = — ( + rac 2 [si n/x Ji-r,(kr) + (-l) n+1 cos /i J_(i_,) (At)] \ , , 

x -- n -^' ) N\ tVS - mc 2 [- sin// J_„(fcr) + (-l) n+1 cos/i J„(Ax)] / 1 J 



with iV a normalization factor Bound state energy I£ again parametrizes different 
self-adjoint extensions and, for given £, determines the parameter [i 



mc 



2 



E {£ 2 - m 2 c A y 



tanu = (-l) n+1 — ^ - ^ • (31) 



By comparing fl3~0|) with (12) one finds a change of the conventional phase shifts 

sin(r77r) . „. 

d = - arctan — - — ( 32 ) 

cos(?77r) + (— l) n tan/x 

When the threshold state £ = mc 2 is present tan \x = and one has a phase shift flip, in 
accord with |T7|]. By the Krein-Friedel formula the contribution of scattering states with 
£ > mc 2 to the integrated density of states is then 

N a {£) - N {£) = -\rj{l - 77) - - arctan - ™ ( ^\ B . ■ (33) 
2 7T cos(r/7r) + (— l) n tan// 

Because of (0) there is again a typical resonance in the relativistic case for < rj < 1/2 
which, however, for 1/2 < r\ < 1 disappear. When r] = 1/2 the resonance is shifted to 
infinity. For negative magnetic field the resonance appears for 1/2 < r\ < 1 and disappears 
for < 7] < 1/2. 

For — £ < —mc 2 the scattering states are given by ^f_g.i(t,r,cj)) = ^^(t, r, 0)| m _»_ m 
supplemented with tan /x — > — tan /i. The former relation is a consequence of /i^lm^-m = 
— /i m for the radial part h m of the Dirac operator and can be verified also by direct 
calculations. The latter ensures the same boundary conditions at the origin as for the 
energy £ scattering states, i.e., the both sets are in the same Hilbert space. Formally, the 
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bound states with the energy ±E and the scattering states with ±/j, are eigenstates of h m , 
but more precise look shows that they belong to different self-adjoint extensions of h m . 

In contrast to the Schrodinger equation it is impossible to find a bound state with 
the homogeneous field regularization and our argument with the decoupling of bound 
states does not work. Whenever the magnetic moment coupling induces an attractive 
V(r) = —g m a/R 2 potential inside the flux tube for the down component then it induces 
a repulsive —V potential for the up component. But it is known that for the bound state 
to exist it has to be at least the attractive potential (EDI). An arbitrary weak attractive 



potential cannot lead to bound states (cf. fll7|). The same situation persists with the 
cylindrical shell regularization. Nevertheless, in the presence of the threshold state the 
total energy of the system with filled Dirac sea changes by infinite positive amount which 
reflects the stability of the system. 

Note that ( p6|) is nothing but radial part of the Klein- Gordon Hamiltonian with the 
Pauli coupling [19| and the relevant results for it are a simple consequence of the above 
calculations. 

4. Discussion, three dimensional case - We have derived the quantum-mechanical cri- 
terion of instability which only involves fundamental parameters of matter and, in 
particular, shows an instability of massless charged particles (cf. PpJ which claims they 



cannot exist in nature as they are completely locally screened in the process of forma- 
tion). We did not find any instability of minimally-coupled relativistic matter without 
the Abelian Chern-Simons term (cf. ]21] which claims that in its presence magnetic field 
is spontaneously generated). 

In 3D the function X(a,g m ) is replaced by X(a, g m )/L, with L the length of the flux 
string. Therefore the above instability may survive up to 3D provided the density of 
states for a long flux ring preserves essential features of the AB potential. However, in 
the case of the 3D Dirac equation any question about the instability due to the symmetry 



of its spectrum is pointless (cf. the suggestion of [22| for 'flux spaghetti' vacuum in the 
spirit of as a mechanism for avoiding the divergence of perturbative QED). With 
the AB potential the pairs of components (x\x 4 )) an d (x 3 ;X 2 ) °f the four-spinor in the 
standard representation combine to the components of two bispinors which satisfy 2D 
Dirac equations with the opposite sign of mass. Because h* m = —h- m for the radial part 
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of the 2D Dirac Hamiltonians (cf. (3) of ||13||), the symmetry of the spectrum (including 
the threshold and bound states) of the 3D Dirac equation with respect to the origin is 



restored, in accord with the general result fIB |. A more general discussion of this and 
related problems will be given elsewhere 



I should like to thank A. Comtet, Y. Georgelin, S. Ouvry, and J. Stern for many useful 
and stimulating discussions. 

References 

[1] We refer to the diamagnetic inequality when the energy of field and matter are 
considered altogether. Sometimes the notion of the diamagnetic and the paramagnetic 
inequality is only used for the energy of matter. See E. Seiler, Gauge Theories as a 
Problem of Contractive Quantum Field Theory and Statistical Mechanics, LNP 159, 
(Springer, Heidelberg, 1982) p. 22 and references therein. 

[2] C. R. Hagen, Phys. Rev. Lett. 64, 2347 (1990); V. V. Semenov, J. Phys. A: Math. 
Gen. 25, L619-L621 (1992); O. Ogurisu, Lett. Math. Phys. 30, 173 (1994). 

[3] P. de Sousa Gerbert and R. Jackiw, Commun. Math. Phys. 124, 229 (1989). 

[4] M. Alford and F. Wilczek, Phys. Rev. Lett. 62, 1071 (1989). 

[5] W. Ehrenberg and R. E. Siday, Proc. Phys. Soc. 62B, 8 (1949); Y. Aharonov and 
D. Bohm, Phys. Rev. 115, 485 (1959); W. C. Henneberger, Phys. Rev. A 22, 1383 
(1980). 

[6] S. N. M. Ruijsenaars, Ann. Phys. 146, 1 (1983). 

[7] M. Abramowitch and I. A. Stegun, Handbook of Mathematical Functions (Dover 
Publ., 1973). 

[8] R. D. Richtmayer, Principles of Advanced Mathematical Physics (Springer, New York, 
1978), vol. 1, ch. 10.15. 

[9] E. C. Marino, B. Schroer, and J. A. Swieca, Nucl. Phys. B200[FS4], 473 (1982). 
There are some misprints and weak points in their calculations but, nevertheless, the 
cited result is valid. 



9 



[10] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Academic 
Press, New York, 1966). 



[11 
[12 

[13 



[14 

[15 
[16 



[17; 



[19 
[20 
[21 
[22 
[23; 
[24 



A. Comtet, Y. Georgelin, and S. Ouvry, J. Phys. A: Math. Gen. 22, 3917 (1989). 
J. S. Faulkner, J. Phys. C: Solid State Phys. 10, 4661 (1977). 

P. de Sousa Gerbert, Phys. Rev. D40, 1346 (1989). Since in the I = -{n + l)-th 
channel the respective order of the Bessel functions for the up/down component is 1 
and when 77 = we have taken v = 1 — r\ instead of v = r\ — 1. 

M. Bordag and S. Voropaev, J. Phys. A: Math. Gen. 26, 7637 (1993). We use e = -|e| 
which leads to I — > — I when compared to their notation. 

Y. Aharonov and A. Casher, Phys. Rev. A19, 2461 (1979). 



We have found that some conclusions of (14| are not valid: their coefficients en, might 



be negative when they claim them to be positive; there is an error in calculating the 
number of bound states. 

C. R. Hagen, Phys. Rev. Lett. 64, 503 (1990). 

B. Thaler, The Dirac equation (Springer, New York, 1992) ch. 5. Unfortunately, he 
does not discuss the bound states. 

I. I. Kogan, Phys. Lett. B262, 83 (1991); J. Stern, Phys. Lett. B265, 119 (1991). 
V. N. Gribov, Nucl. Phys. B206, 103 (1982). 

Y. Hosotani, Phys. Lett. B319, 332-318 (1993); UMN-TH-1 238/94 (|hep-th/940209q ). 



A. S. Goldhaber, H. Li, and R. R. Parwani, ITP-SB-92-40 (|hep-th/9305007| ). 



H. B. Nielsen and P. Olesen, Nucl. Phys. B160, 380 (1979). 
A. Comtet, A. Moroz, and S. Ouvry, IPNO/TH 94-30. 



10 



